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ABSTRACT 

Evidences of non-thermal X-ray emission and TeV gamma-rays from the su- 
pernova remnants (SNRs) has strengthened the hj^othesis that primary Galactic 
cosmic-ray electrons are accelerated in SNRs. High energy electrons lose en- 
ergy via synchrotron and inverse Compton processes during propagation in the 
Galaxy. Due to these radiative losses, TeV electrons liberated from SNRs at dis- 
tances larger than ~1 kpc, or times older than ~ 10^ yr, cannot reach the solar 
system. We investigated the cosmic-ray electron spectrum observed in the so- 
lar system using an analytical method, and considered several candidate sources 
among nearby SNRs which may contribute to the high energy electron flux. Es- 
pecially, we discuss the effects for the release time from SNRs after the explosion, 
as well as the deviation of a source spectrum from a simple power-law. From this 
calculation, we found that some nearby sources such as the Vela, Cygnus Loop, 
or Monogem could leave unique signatures in the form of identifiable structure in 
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the energy spectrum of TeV electrons and show anisotropics towards the sources, 
depending on when the electrons are liberated from the remnant. This suggests 
that, in addition to providing information on the mechanisms of acceleration and 
propagation of cosmic-rays, specific cosmic-ray sources can be identified through 
the precise electron observation in the TeV region. 

Subject headings: cosmic-rays electrons — supernova remnants — propagation 
— acceleration 



1. Introduction 

Radio observations have indicated that supernova remnants (SNRs) are the most likely 
sources of cosmic-ray electrons in the energy region below ~ 10 GeV. Evidence for non- 
thermal X-ray emission from the supernova remnant SN1006 discovered with the ASCA 
satellite (Koyama et al. 1995) strongly supports the hypothesis that Galactic cosmic- 
ray electrons in the TeV region originate in supernovae. In this case, TeV gamma-rays 
should be produced via the inverse Compton process between accelerated electrons and the 
cosmic microwave background (CMB) radiation, and indeed TeV gamma-rays were detected 
by the CANGAROO experiment (Tanimori et al. 1998). Additional evidence for X-ray 
synchrotron emission is provided by observations of several other SNRs, such as RX J1713.7- 
3946 (Koyama et al. 1997), Gas A (Allen et al. 1997), 1G443 (Keohane et al. 1997), 
G374.3-0.5 (Slane et al. 1999), and RX J0852.0-4622 (Slane et al. 2001). In SNR RX 
J1713. 7-3946, evidence for the acceleration of cosmic-ray protons is reported by Enomoto 
et al. (2002) from TeV gamma-ray observations. However, there are also arguments that 
the TeV gamma-rays observed cannot be interpreted as hadronic in origin (Reimer & Pohl 
2002; Butt et al. 2002). 

To clarify the origins of cosmic rays and their propagation mechanisms within the 
Galaxy, electrons provide an ideal probe, complementary to the nucleonic component, due 
to their low mass and leptonic nature. High energy electrons lose energy primarily via syn- 
chrotron and inverse Gompton processes during propagation in the Galaxy. These processes, 
combined with the absence of hadronic interactions, simplify modeling of the propagation 
of electrons compared with other cosmic-ray components such as nucleons. Since the energy 
loss rate is almost proportional to -E^, higher energy electrons lose energy more rapidly. TeV 
electrons lose most of their energy on a time scale of ~ lO^yr, and their propagation distances 
are therefore hmited to several hundred pc. 

Measurements of cosmic-ray electrons in the TeV region have been successfully per- 
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formed only with emulsion chamber detectors (Nishimura et al. 1980; Kobayashi et al. 
1999). The observed energy spectra extend without cut-off up to ~ 2TeV. This means that 
the observed TeV electrons must have been accelerated in SNRs at distances within several 
hundred pc, and at times within ~ 10^ yr ago. 

Several models have been proposed to describe cosmic-ray propagation in the Galaxy, 
motivated by attempts to fit existing data on the heavy primary composition and spectra 
of different species assuming the same simple power-law source spectrum. Among these 
models, it is believed that the diffusion model provides the most realistic description of 
the propagation. The solution for the electron density in the diffusion equation is derived 
by several authors (e.g. Jokipii & Meyer 1968; Ginzburg & Ptusk;in 1976; Nishimura et 
al. 1979; Berezinskii et al. 1990 and other references therein), under different boundary 
conditions for the Galactic halo. 

In a study of the propagation of cosmic-ray electrons, Shen (1970) first pointed out 
that a continuous source distribution model is not valid in the energy region above lOOGeV, 
because the electron spectrum in that energy range depends upon the age and distance of a 
few local sources. He showed that high energy electrons over lOOGeV probably come from 
the Vela X pulsar, and predicted a cutoff at ~2 TeV in the energy spectrum. Although 
various parameters for SNRs and diffusion coefficients were not known very well at that 
time, the concepts proposed by the author has been accepted in later work. 

Cowsik & Lee (1979) have investigated the electron spectrum in terms of contributions 
from discrete sources such as SNRs distributed over the Galaxy. Their calculations apply 

in the the case of sources continuously active in time. They suggested that such sources 
must be located within a few hundred pc of the solar system to prevent radiative energy 
losses from inducing a premature cut-off in the energy spectrum. They took the diffusion 
coefficient D = 10^^ cm^s~^, independent of energy, which is much smaller than recently 
accepted values in the energy region over 100 GeV (see section 2.2). They concluded that it 
is very unlikely that SNRs are the only sources of cosmic-ray electrons in the energy range 
IGeV — ITeV, but this result seems due to the inappropriate value of the diffusion coefficient 
assumed by themselves. 

Nishimura et al. (1979) calculated the propagation of electrons from SNRs, using a 
solution of the diffusion equation with a disk-shaped boundary condition. They showed 
that the TeV electron spectrum would deviate from power-law behavior due to fluctuations 
caused by the small number of sources capable of contributing to the observed flux in the 
TeV region. 

To explain the features of both energy spectrum and charge composition of electrons. 
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Aharonian et al. (1995) suggested an approach that separates the contribution of one nearby 
source, with age lO^yr and distance lOOpc, from the contributions of other Galactic sources at 
distances over Ikpc. Using this approach, Atoyan et al. (1995) calculated the electron energy 
spectrum in analytical form. They considered the energy-dependent diffusive propagation 
of electrons, assuming the diffusion coefficient D = 7 x lO^^cm^s"^ at 10 GeV, increasing 
with energy as E^-^. They showed that the observed energy spectrum from sub-GeV to TeV 
energies can be explained by this model. 

Pohl & Esposito (1998) has pointed out that, if cosmic-ray electrons arc accelerated 
in SNRs, the observed electron spectrum should be strongly time-dependent above 30 GcV, 
due to Poisson ffuctuations in the number of SNRs within a given volume and time interval. 
They assumed the diffusion coefficient D = 4xlO^''(E/lGeV)°''cm^s~^, and calculated 400 
cases of spectra using a random distribution of SNRs in space and time. They argued that 
in their model, an injection spectral index of 2.0 cannot be rejected for the observed local 
electron flux at high energy, with a probabihty of 5% over all their calculated spectra. 

Erlykin & Wolfendale (2002) considered the model for cosmic-ray electron production 
in the supernova explosion and variants for the subsequent propagation. In their model, 
electrons are accelerated throughout the SNR life time of several 10^ yr, and after the cease 
of the acceleration process electrons leave the remnant. They simulated the electron spectra 
by taking the random distribution of SN in the Galaxy and letting the electrons diffuse to 
the solar system. 

As pointed out in these works, at energy higher than ~1 TeV, the propagation lifetime 
of electrons is so short that only a few cosmic-ray sources can contribute to the observed 
ffux. Thus one expects ffuctuations in the local electron energy spectrum reffecting the 
contributions of these small number of discrete sources. To analyze ffuctuations in the 
spectrum, we use the method of separating the contributions of distant and nearby sources 
to the total ffux of high energy electrons. Although this approach was discussed by earlier 
authors (e.g. Atoyan et al. 1995 and references therein), the main difference between the 
present and earlier works is that in this paper we took only known, observed local SNRs 
as nearby sources. Using information on the ages and distances of observed SNRs in the 
neighborhood of the solar system, we could determine which sources contribute electrons 
most efficiently in a given energy range, and thus correlate individual sources with features 
of the electron spectrum. In addition, one would expect the high energy electron ffux to be 
anisotropic if the most signiffcant sources are nearby (Shen & Mao 1971; Ptuskin & Ormes 
1995). This means that we can identify cosmic-ray electron sources from the analysis of the 
local electron spectrum together with anisotropics in arrival direction. 

For this purpose, we calculated the electron energy spectra with the diffusion model in 
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a semi-analytic approach, separating the contributions of distant and nearby sources. In this 
calculation, we particularly took into account of the effects of the delay of electron liberation 
from SNRs after the explosion, as well as the deviation from a single power-law spectrum 
of electrons in SNRs. As a result, we indicate that nearby SNRs such as the Vela, Cygnus 
Loop, or Monogem are the most likely candidates for the sources of cosmic-ray electrons in 
TeV region, depending on when the electrons are released from SNRs. If these sources are 
proven to be the main contributors of high energy electrons, we can precisely analyze the 
propagation of TeV electrons by comparing predicted results with observed electron spectra, 
and the degree of anisotropy in the direction of the sources. Such analyses will yield more 
detailed information about the acceleration and propagation of cosmic rays. 

2. Acceleration in SNRs and Propagation in the Galcixy 

2.1. Accelerated Electrons in SNRs and their Liberation Time 

According to shock- acceleration models, the maximum energy of accelerated electrons 
is limited, by the SNR age, a free escape from the shock region of MHD turbulence, or 
synchrotron losses (Sturner et al. 1997; Reynolds 1996). Analysis of the observed radio 
and X-ray spectra also indicates that the typical electron spectrum produced in a remnant 
is a power-law with a cut-off of ~10 — lOOTeV (Reynolds & Keohane 1999; Hendrick 2001). 
Therefore, we approximated an electron injection spectrum of a power-law with an exponen- 
tial cut-off of the form of exp(—E/Ec), as taken by these authors. 

As discussed in Erlykin & Wolfcndale (2002), it is also important to take account of 
when the accelerated electrons arc liberated from the remnant after the explosion. From 
shock-acceleration theory, it is suggested that the electrons are liberated from SNRs in the 
termination of the shock, i.e. when the shock velocity has dropped to the mean ISM Alfven 
velocity and the remnant has been reduced in the ISM (e.g. Dorfi 2000). This time scale 
is ~ 10^ yr. However, there are no clear observational evidences for the release time of 
the TeV electrons from SNRs. In the scenario of the SNR origin of cosmic-ray electrons, 
the accelerated TeV electrons should be liberated from the remnant within ~10^yr, or the 
acceleration process should operate till the last phase of SNRs. This is because the life time 
of TeV electrons in the remnant is estimated to be less than a few times of 10*^yr due to 
synchrotron loss inside SNRs with a magnetic field of 10//G. Therefore, we evaluated the 
cosmic-ray electron spectra, changing the electron release time r from to 1 x 10^ yr. 
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2.2. Properties of the Propagation of High Energy Electrons 

High energy electrons above 10 GeV lose their energy mainly via the synchrotron and 
inverse Compton processes while propagating through the Galaxy. The energy loss rate is 
given by 

with 

Here, E is the electron energy, m is the mass of electron, c is the speed of light, B is the 
magnetic field strength in the Galaxy, is the energy density of interstellar photons, and 
(T is the cross-section for Thomson scattering. This formula (2) is valid when the inverse 

Compton process is well approximated by the Thomson scattering cross-section, but we need 
corrections from the Klein-Nishina formula for high energy electrons as described below. 

Typically quoted values for the interstellar magnetic field are estimated from measure- 
ments of the Zeeman splitting, the Faraday rotation, or the radio synchrotron emission. The 
estimates from Zeeman splitting and Faraday rotation measurements refer only to the line 
of sight component of the magnetic field. They are also biased toward cold neutral and 
warm ionized regions, respectively. Using the radio synchrotron emission from relativistic 
electrons, the local magnetic field strength of S_l~5 /iG is obtained as a more adequate value 
(Ferriere 2001). Here, B± means the magnetic field perpendicular to the electron velocity, 
that is = 2S^/3. In this calculation, we adopted B± — 5//G for our Galaxy. 



Figure 1 



Since the Thomson scattering cross-section approximation is inadequate for high energy 
electrons, we evaluate the energy loss coefficient b using the Compton scattering cross-section 
in the following. The interstellar radiation field in the Galaxy is dominated by three com- 
ponents: 2.7K cosmic microwave background (CMB), re-emitted radiation from dust grains, 
and stellar radiation. The energy densities of photons are 0.26 eV/cm^ for CMB, 0.20 eV/cm^ 
for re-emitted radiation from dust grains, and 0.45 eV/cm^ for stellar radiation, respectively 
(Mathis et al. 1983). For inverse Compton scattering, we calculated the Klein-Nishina 
cross-section accurately using the formula by Blumenthal & Gould (1970). The resulting 
energy loss coefficient b decreases gradually with energy as shown in Figure 1, where the ra- 
diation field is assumed to be isotropic. For electrons in the TeV region, energy loss through 
interactions with the CMB dominates over the other two components. 
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In a diffusive propagation model, the diffusion coefficient determines the travel distance 
of electrons in a given time. The diffusion coefficient D has been estimated by using the 
observed ratios of secondary to primary nuclei (B/C) by HEAO-C (Engelmann et al. 1990) 
and Voyager (Lukasiak et al. 1994), and is given by 

L> = 2 X 102^(£;/5GeV)^(cm7s), (3) 

where 

r {E < 5GeV) 
~\0.6 (E>5GeV). 

The formula (3) agrees well with the experimental data up to tens of GeV. Beyond lOOGeV, 
there have been scarcely data on B/C until quite recently, and it was not clear whether the 
formula (3) is applicable to the higher energy region or not. 

There are arguments that simple extrapolation of DocE^ with 5 = 0.6 gives too large 
diffusion coefficient in 1 — 100 TeV region to interpret the observed anisotropy of primary 
cosmic-rays less than 10~^ (Ambrosio et al. 2003 and other references therein). Thus, it has 
been argued that the increase of the diffusion coefficient with energy should become slower 
than 5 — 0.6 at somewhere beyond tens of GeV, and also that 6 — 0.3 is plausible from 
a Kolmogorov-type spectrum of turbulence in the interstellar medium (e.g. Gaisser 1990, 
2000; Jones et al. 2001). In fact, a recent work of emulsion chamber experiments with 
long balloon exposures has shown the data of B/C and (secondary species of Fe)/Fe in TeV 
region, which indicate that the diffusion coefficient is consistent with the change of S from 
0.6 to 0.3 between tens of GeV and ITeV (Furukawa et al. 2003). 

Hence we take D = (2 — 5) x 10^^(£'/TeV)°-^(cm^s"^) as a plausible diffusion coefficient 
in TeV region. The numerical values of 2 and 5 are derived from the assumption that the 
formula (3) is valid up to 50 GeV or 1 TeV, respectively. With this assumption, the source 
spectral index 7 is chosen to maintain 7 + 5 = 2.7. 

Since high energy electrons lose energy by synchrotron and inverse Compton processes 
at the rate dE/dt = —hE"^, electrons lose almost all of their energy E after time T = 1/bE. 
Therefore, electrons observed with energy E must have been accelerated within T = 1/bE 
from the present. Hence, the lifetime T becomes progressively shorter with increasing energy. 
Assuming B± = 5/iG and taking the Klein- Nishina formula for Compton process, the lifetime 
is T = 1/bE — 2.5 x 10^(yr)/£'(TeV). As an approximate treatment, electrons can diffuse a 
distance of = (2Dr)V2 during this time; i.e. 0.6 - 0.9kpc for D = (2 - 5) x 1029(cm2s-i) 
at 1 TeV. More precise estimate for R is given in the following section. 
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2.3. Propagation of Cosmic Rays from a Single Source 

In the diffusion model for tlie propagation of electrons in the Galaxy, the electron density 
Ne is given by the equation 

^ - VpVTVe) - ^(bE'N,) = Q(E,r,z,t), (4) 

where Q is the electron source strength, and r is the distance to sources from the solar 
system. 

We shall take into account the fact that the Galactic disk, where the sources are dis- 
tributed, is surrounded by a halo in which the cosmic rays are confined for a long time before 
they escape into intcrgalactic space. We assume that cosmic-ray electrons are bounded by 
two parallel planes at z = ±h with the median plane being occupied by the Galactic disk. 
The halo thickness, h, is estimated to be 2.8lQ;gkpc from ^°Be observations from the HET 
aboard the Voyager 1 and 2 spacecraft (Lukasiak et al. 1994), which is also consistent with 
the observed Galactic radio emission structure (Beuermann et al. 1985). Thus we need 
to consider solutions with these boundary conditions for further analysis. Since the solar 
system is located at 2; = 15 pc, quite near the median plane of the disk, and the sources are 
uniformly distributed in layer much thinner than the confinement domain, we put z = 0. 
Assuming burst-like injection of electrons after the supernova explosion, the source term can 
be represented by 

QiE,r,t)=Q{E)S{r)5{t). (5) 
The density of electrons, N^, from a point source with injection spectrum 

Q(E) = QoE-^eM-E/E,) (6) 

at a distance r and time t after the release of electrons from the remnant, is derived from 
the diffusion equation using the Fourier transform, assuming the diffusion coefficient has the 
form D — DqE^. Taking the boundary condition A^e = at the boundary of the Galactic 
halo z — ±/i, the general solution of the equation (4) in cylindrical coordinates is given by 

°° O 

N,{E,r,t) = Yl j^exp(-L'iA;2-rV(4L'i))cos(A;„^)(l-6£;t)^-2£;-^exp(-£;/(£;e(l-&£^i))) 

n=0 ^ 



where 



and 



(7) 



Di^ {l/b)DE-^dE^ 
Je 



6(1 - 6)E^-^ 
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2.4. Distant Components and Nearby Components 



Integrating Ne{E,r,t) of the solution (7) with r from to oo, we derive 



N^{E,t) = / N^{E,r,t)27rrdr 



= J2 l^^M-Dikl){l - bEty-'E-^eM-E/iE^il - bEt))). (8) 



In this calculation, we approximated the integration on r to range from to oo. Since 
electrons cannot propagate to larger distances than the halo thickness of h in the radial 
direction, we may ignore the effects of the lateral distribution of the sources on the Galactic 
disk. In integrating Ne{E,t) of the formula (8) on t from to 1/bE, we need the numerical 
integration. In the case of << Ec, however, we can perform the integration with an 
analytical way, since the exponential cut-off term can be ignored. The solution is given by 
a Confluent Hyper-geometric Function for each term in the series. 



and / is the supernova explosion rate in the Galaxy per unit area and unit time. This 
analytic form is convenient to sum up many terms in the series. Higher energy electrons 
cannot reach the halo boundary due to their larger energy loss rate, and distribute around the 
Galactic disk plane. Therefore, we need to include the higher-order Fourier components to 
estimate the flux of higher energy electrons. In fact, we need to add up to 10 terms at lOGeV 
while up to 2 X 10^ terms at lOTeV to get a precision of a few percent in the convergence 
of the series solution (9). In order to simplify the calculation for the high energy region, we 
can use a three dimensional solution without boundary conditions instead of using the two 
dimensional solution (9), which is given by (e.g. Berezinskii et al. 1990 and other references 
therein) 





where iFi is a Confluent Hyper-geometric Function and 




Ne{E,r,t) 



1 




E 



)(1 - bEt^eM-E/iE^il - bEt))), (10) 



(47rDi)3/2 



1-bEt 



where 
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The accuracy obtained by this approximation is discussed in Appendix A. 

The flux of electrons from sources continuously and uniformly distributed in the Galactic 
disk is derived from the integration of electron density with respect to r and t in the following, 

J{E) = ^N^{E). (11) 

In the low energy region, below 1 TeV, many sources contribute to the observed electron flux. 
As the energy of the observed electrons increases, the number of electron sources decreases 
and only nearby sources can contribute to the electron flux. As described in section 2.2, 
young cosmic-ray electron sources which are nearby create discrete effects in the TeV region, 
such as structure in the electron spectrum and the anisotropy favoring the source direction. 
Thus we can deflne the contribution from distant or relatively old sources by subtracting the 
flux of nearby young sources, eliminating the effects of fluctuations due to the small number 
of nearby sources, as follows: 

UE) = JiE)-UE), (12) 

where Jd and Jn show the flux of electrons from distant (in space-time) sources and nearby 
sources, respectively. Here, we calculated the contribution of nearby sources distributed 
uniformly in space and time as 

UE) = —J^ dtj^ drfN,{E,r,t)27rr. (13) 

Appropriate values of Rq and Tq are discussed in section 4.2. 



3. SNRs as the Cosmic- Ray Electron Sources 

As described in section 1, recent X-ray and TeV gamma-ray observations indicate that 
high energy electrons are accelerated in SNRs. Here, we assumed that SNRs are cosmic-ray 
electron sources, and that electrons are hberated burst-likely from SNRs at the time of r = 0, 
5x10^ yr, 1x10^ yr, 5x10^ yr, and 1x10^ yr after the explosion. As described in section 4, 
the case of the continuous release is well approximated by the burst-like release with a mean 
value of the liberation time. 

Prom the observations of SN1006 (Reynolds 1996; Koyama et al. 1995; Tanimori et 
al. 1998), the output energy in electrons above 1 GeV is estimated to be ~ 1 x 10^^ erg. 
We can also estimate the net output energy for electrons from supernovae using the energy 
density of cosmic rays and the supernova rate of l/30yr in the Galaxy (Berezinskii ct al. 
1990; Gaisser 2000 and other references therein). The output energy of electrons above 1 
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GeV estimated in this way is about 10^^ erg/SN, and is consistent with the estimate from 
SN1006. Therefore, we took the output energy from a supernova in electrons above 1 GeV 
is W ^IxW^ erg. 

To analyze contributions due to discrete cosmic-ray electron sources, we list in Table 1 
all known SNRs which are located within 1 kpc of the solar system and with age less than 
4 X 10^ yr. For each source, the maximum energy E'max of electrons reaching the solar system 
is given by -Emax = ^/{bT), with the energy loss coefficient b shown in figure 1 and the age 
T of each SNR, which are determined independent of the diffusion coefficient. 



Table 1 



Distance to SNRs is an important parameter in evaluating the high energy electron flux, 
and we comment on this parameter of some SNRs listed in Table 1. 

So far, the canonical distance to the Vela SNR has been taken to be 500pc, a value which 
was derived from the analysis of its angular diameter in comparison with the Cygnus Loop 
and IC443 (Milne 1968), and pulsar dispersion determination (Taylor & Cordes 1993). 
However, recent parallax measurements clearly indicate that the distance of 500pc is too 
large. Cha et al. (1999) obtained high resolution Ca-II absorption fine toward 68 OB stars in 
the direction of the Vela SNR. The distances to these stars were determined by trigonometric 
parallax measurements with the Hipparcos satellite and spectroscopic parallaxes based upon 
photometric colors and spectral types. The distance to the Vela SNR is constrained to be 
250 ± 30pc due to the presence of the Dopplcr spread Ca-ll absorption line attributable to 
the remnant along some lines of sight. Caraveo ct al. (2001) also applied high- resolution 
astrometry to the Vela pulsar (PSR B0833-45) V ~ 23.6 optical counterpart. Using Hubble 
Space Telescope observations, they obtained the flrst optical measurement of the annual 
parallax of the Vela pulsar, yielding a distance of 2941^0 pc. Therefore, we calculate the 
electron flux adopting a distance of 300 pc to the Vela SNR. 

Previously, a distance of 770 pc to the Cygnus Loop was often quoted (Minkowski 1958). 
Recently Blair et al. (1999) suggested that the distance is 440;!;}qq pc, based on Hubble Space 
Telescope observations of a fllament in the remnant. Thus we took the distance to be 440 
pc, considerably smaller than the previously quoted distance of 770 pc. 

For the recently discovered SNR RX J0852.0-4622, the estimated age of ~ 680yr and 
distance of ~ 200pc was based on thermal X-ray emission and ^^Ti gamma-ray detection 
(Aschenbach 1998; lyudin et al. 1998). However, ASCA observations reveal that the X- 
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ray emission is non-thermal, and re-analysis of the COMPTEL data indicates that the Ti 
gamma-ray detection is only significant at the 2 — Act level (Slane et al. 2001 and references 
therein). Thus it is suggested that RX J0852. 0-4622 is at a larger distance, 1 — 2kpc (Slane 
et al. 2001), so we do not include it in Table 1. 

4. Results 

4.1. Electrons from the Known Necirby SNRs 

Figure 2 shows contours of the expected electron flux at 3 TeV (with flux values scaled 
by ) as a. function of age and distance of the SNR. As described in the preceding section, 
we assumed that the output energy of electrons over 1 GeV is FT = 1 x 10^® erg/SN, and 
the injection spectrum is a power-law function with an exponential cut-off of the form of 
e-x-p{—E/Ec). The diffusion coefficient has the form D — Do{E /TeVY cm^s~^ where 5 is 
0.3 in the TeV region and the spectral index 7 is chosen to maintain ^ + 5 — 2.7. As shown 
in Figure 2, the electron flux is strongly dependent on source age and distance. In the case 
of the prompt release of electrons after the explosion, the flux from the Vela SNR is the 
largest among the known SNRs listed in Table 1. The flux value is quite sensitive to the 
change of distance to Vela from 500pc to 300pc, since the solution for electron density yields 
a Gaussian distribution function of r as shown in formula (7) or (10). The flux of electrons 
at a distance of 300pc is two orders of magnitude larger than at 500pc. 



Figure 2 



4.2. Calculated Electron Energy Spectrum Compared with the Observed 

Spectrum 

We separately calculated the contributions to the electron energy spectrum from nearby 
and distant sources. We selected SNRs in the neighborhood of the solar system located within 
a distance of Ikpc and an electron release time within lO^yr in the past. The boundary 
condition in this domain almost corresponds to the diffusion path length and life time of 
electrons around 1 TeV. We calculated the electron flux from nearby sources using the 
distance and ages of selected SNRs in Table 1 with formula (10). Justification of taking 
Rq = Ikpc and Tq = lO^yr is discussed in Appendix B. The SNRs in the nearby region are 
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thus SN185, S147, HB 21, G65. 3+5.7, Cygnus Loop, Vela, and Monogem. In this calculation, 
we assumed that supernovae occur uniformly on the Galactic disk at the rate of 1/30 yr, 
and took the halo thickness tohe h — 3kpc. 

Figure 3 shows the calculated energy spectra of electrons without a cut-off of the injec- 
tion spectrum in the case of the prompt release after the explosion {t — 0), compared to 

the observed data (Golden et al. 1984; Tang 1984; Golden et al. 1994; Kobayashi et al. 
1999; Boezio et al. 2000; Alcaraz et al. 2000; DuVernois 2001; Torii et al. 2001). Here, 
we illustrate the cases of the diffusion coefficient of -D = Dq{E /TeV)^''^ with Dq = 2 x 10^^ 
crn^s"^ and Dq = 5x 10^^ cm^s~^ in TcV region, with D = 2x 10^^(£'/5GcV)°-^ cm^s"-^ upto 
50 GeV or 1 TeV as given by the formula (3) in section 2.2. In this figure, we also plotted the 
flux of low energy interstellar electrons estimated from the Galactic radio data (Rockstroh 
& Webber 1978). Cosmic-ray electrons have been observed by a variety of instruments, but 
only emulsion chamber data provide observations of electrons above 1 TeV. Since the flux 
of TeV electrons is low, detectors must have large acceptance Sfl and high proton rejection 
power (~ 10^). The emulsion chamber satisfies these requirements (Nishimura et al. 1980; 
Kobayashi et al. 1999). 

Our propagation model, which consists of separately calculated distant and nearby 
components, is consistent with the observed data for the local primary electron spectrum in 
the energy range from lOGeV to 2TeV. Our model also predicts that nearby SNRs insert 
unique, identifiable structures in the electron spectrum from ITeV to lOTcV. As shown in 
Figure 3, the absolute fiux and spectral shape change with the diffusion coefficient. On the 
other hand, the maximum energy of each SNR is the same, independent of the diffusion 
coefficient value, because it is determined by the age of the SNR. 



Figure 3 



We also calculated how the source spectral shape and release time of electrons from 
the remnant affect the cosmic-ray electron spectrum. We illustrate some examples of these 
calculations by taking Dq = 2 x 10^^ cm^s^^. Figure 4 shows the calculated energy spectra 
with a cut-off of =10 TeV, 20 TeV, and oo in the injected electron spectrum, assuming 
the prompt release after the explosion. We can find that these spectra are similar with each 
other, independent of the cut-off energies. Figure 5 shows the calculated energy spectra with 
a cut-off of Ec — 20TeV, in which electrons are released burst-hkely after the explosion in 
the release time of r = 5 x lO^yr, 1 x lO^yr, 5 x lO^yr, and 1 x lO^yr, respectively. As shown 
in this figure, it is to be noted that the spectrum for r = 5xl0^yr is almost the same with 
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that of the prompt release after the explosion. However, the delay of the release time from 
SNRs have a large impact on the flux in the TeV region for T>lxlOVr. 



Figure 4 



Figure 5 



We also checked the case of the continuous release. Figure 6 shows that these spectra 
are well represented by that of the burst-like release with a mean value of the continuous 
release time. 



Figure 6 



5. Summary and Discussion 

We calculated the energy spectrum of cosmic-ray electrons, separating the contributions 
of sources nearby in space and time (those within distances of Ikpc and times of 1 x lO^yr) 
and Galactic sources located outside this domain. 

There are only 9 SNRs within a distance of Ikpc and an age of 4 x lO^yr, as shown in 
Table 1, which is much smaller than the expected number of ~ 60 assuming a supernova 
rate of l/30yr and a Galactic disk radius of 15kpc. This is due to selection biases in radio 
observations, since surface brightnesses in radio observations decrease with age. Thus, also 
coupled with the reason of the arguments in Appendix B, we define a domain of R <lkpc 
and T < 1 X lO^yr to evaluate the electron flux from nearby sources. There are 7 known 
SNRs within this region listed in Table 1: Vela, Cygnus Loop, Monogem, G65.3-I-5.7, HB 
21, S147, and SN185. On the other hand, the expected number of SNRs in this smaller 
domain is ~ 15. The difference could be partly due to statistical fluctuations in the small 
number of sources, and it could also be due to undetected SNRs in this domain. In fact, 
the surface brightness estimated from an age of 1 x lO^yr, assuming adiabatic phase, is still 
fainter than the typical detection limit of surface brightness in studies of the distribution 
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of Galactic SNRs (Kodaira 1974; Leahy & Xinji 1989). Although this indicates that the 
electron flux calculated from known nearby SNRs may give a lower limit of the flux, the 
contribution of undetected SNRs is expected to be relatively small in the case of the prompt 
release of electrons after the explosion. Since these undetected SNRs should have lower 
surface brightness and be older than detected SNRs, they could scarcely contribute to the 
electron intensity in the TeV region as indicated in figure 2. Some quantitative arguments 
on this point are discussed in Appendix B. 

Our model predicts that nearby SNRs such as the Vela, Cygnus Loop, or Monogem, 
present unique, identifiable structures in the electron spectrum from ITeV to lOTcV. As 
the diffusion coefficient increases, electrons propagate larger distances with longer mean 
path length, resulting the electron density becomes smaller. This causes a flatter spectral 
shape and smaller peak flux, as shown in flgures 3. The release time r after the explosion 
determines which SNRs contribute for cosmic-ray electrons in the TeV region. For r — 
and T = 5 X lO^yr, the energy spectra are similar with each other, and the Vela SNR is the 
most dominant source in the TeV region. For r = 1 x lO^yr and r = 5 x lO^yr, the Cygnus 
Loop and the Monogem SNR are dominant, and for r = 1 x lO^yr there are no dominant 
known sources in the TeV region. We can see that some combinations of parameters (e.g. 
D = 5 X 10^^(£^/TeV)°'^ cm^s~^, r = 0, and = oo) are already ruled out even by the 
limited existing data. 

Besides such primary electrons accelerated in SNRs, the observed cosmic-ray electron 
flux includes secondary electrons, produced mainly by interactions of cosmic-ray protons 
and nuclei with interstellar gas. The secondaries are mostly decay products of charged pions 
produced in interactions, i.e. tt^ — > //^ — > e^. As positrons are produced by this process, 
we can estimate the intensity of secondaries relative to primaries to be ~ 10% at 10 GeV 
from the observed ratio of e+ to e~ + e+ ~ 5% at 10 GeV (Barwick et al. 1997; Boezio 
et al. 2000). The spectral index 7 = 2.7 for secondaries is the same as the index of the 
parent cosmic-ray protons and nuclei. This index is larger than that of primary electrons 
at the source, which is ~2.4 in the TeV region. Since the intensity relative to primary 
electrons decreases with increasing energy as (£'/10GcV)~^'^~^''' = (E/lOGcV) the 
relative intensity of secondaries becomes ~2.5% at 1 TeV. This suggests that the intensity 
of the secondary electrons is negligible compared to primary electrons in the energy range 
of interest here. 

Next, we estimate the degree of anisotropics from the speciflc sources in TeV region. 
Let Ni[E,ri,ti) be the contribution to the local density of cosmic-ray electrons at energy E 
from a source located at distance and of age tj. The anisotropy parameter Aj due to the 
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density gradient of electrons is given by 




3£> VNj 



2cti 



3ri 



(14) 



for an individual source, where /max and /min are the maximum and minimum electron 
intensity in all directions (Shen & Mao 1971; Ptuskin & Ormes 1995). For example, the 
anisotropy of electrons for Vela with r = is estimated to be 13%, where n = SOOpc and 
ti — 1.1 X lO^yr. The anisotropy becomes larger if we include the effect of the delay of the 
release time from SNRs. 

In this paper, it is demonstrated that measurements of the energy spectrum of electrons 
in the TeV region are crucial to detect the unique effects of nearby sources as described above. 
If we observe pronounced features in the shape of the spectrum, together with anisotropy 
towards the nearby SNRs, we would confirm the nearby SNRs as the main contributors to 
electrons in the TeV region. We can also perform a more detailed analysis for the release 
time of electrons, cut-off energy, output energy, diffusion coefficient, and so on. 

High energy cosmic-ray electrons are the most powerful probe to help us identify the 
origin of cosmic rays. However, at present, only emulsion chamber data are available in the 
energy region 1 ~ 2TeV. To make significant additional observations on electrons in this 
region, not only balloon flight experiments with long exposures but also new experimen- 
tal programs such as the CALET (CALorimetric Electron Telescope) on the International 
Space Station (Torii et al. 2002) are to be promising to observe the electron spectrum 
and anisotropy in TeV region with high statistical accuracy. Such direct observations will 
reveal the origin of cosmic-ray electrons in this energy region, and also bring us important 
information on the sources, acceleration, and propagation of cosmic-ray electrons. 

We are grateful to Prof. R.J.Wilkes for his careful reading and useful comments of the 
manuscript. 

A. Relation between the three dimensional solution of the diffusion equation 
without boundaries, and the two dimensional solution with boundaries 

In section 2.4, we gave the solution (9) of the two dimensional diffusion equation for a 
source of a power-law electron spectrum with boundaries defined by parallel planes. How- 
ever, since it has the form of a series, it is tedious for numerical evaluation and hard to 
understand the physical meaning of each term. We show that we can simplify the situation 
at high energies by substituting a three dimensional solution without boundaries for the two 
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dimensional solution with boundaries. In this appendix, to see a general behavior of the 
solution, we treat a source of a power-law spectrum without an exponential cut-off. The 
arguments are vahd for E « E^, that is £; < ~lTeV for E^ ^ 10 - 20TeV. 

Electrons propagate an average distance oi R — y/2DT during their lifetime of T = 
1/bE. If the propagation distance is smaller than the galactic halo thickness h, a three 
dimensional solution without boundaries is applicable. Assuming that a source is located on 
the plane of the Galactic disk, the solution of the diffusion equation (4) is given by 

where 

■'E/(l-bEt) 



Di= {l/b)DE-^dE. 
Je 



Here, assuming that the diffusion coefficient has the form D — DqE^ and a source spectrum 
of the form Q{E) — QqE^'^ , we have (e.g. Ginzburg & Ptuskin 1976 and other references 
therein) 

= (iJ^""''^^'"''^^^ " hEtr~'E-\ (A2) 

where 

_ L>o(l - (1 - hEtf-') 
' 6(1 - 5)E^-^ 

Integrating the solution Ne{E, r, t) (A2) on r from to cxd and on t from to 1/bE, the 
electron flux J3 (from sources assumed to be continuously and uniformly distributed in the 
Galactic disk) is given by 

where B is the Beta function. 

We compared this three dimensional solution of J3 with the two dimensional solution 
of J with boundaries described in formula (9). Table 2 shows ratios of J3 to J for halo 
thicknesses of Ikpc to 5kpc and electron energies of IGeV to lOTeV, taking the case of 
a diffusion coefficient of D = 2x 10^^(i?/TeV)°'^cm^s"^ as an example in the TeV energy 
region with formula (3) below 50GeV. For larger halo thicknesses or electron energies, the 
differences between J3 and J become smaller. For the case where halo thickness h = 3kpc, 
it is shown that the two solutions above lOGeV region agree with each other within an error 
of 1%. 
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Table 2 



B. Evaluation of Cosmic-ray Electron Flux from Distant Sources 

We define the total fiux of electrons as J{E), the fiux from distant source at large 
distances or relatively old ages as Jd, and the flux from nearby sources as Jn. 

To avoid statistical fluctuations in the flux due to the small number of nearby sources, 
we evaluate Jd as shown in section 2.4 in the following, 

ME)^J{E)-UE), (Bl) 

where 

Jn = ^y dtj drfN^{E,r,t)27Tr (B2) 

is formula (13) defined in section2.4. 
Then, we have 



poo poo j-lQ j-tto 

ME) = ^{y dtj drfN,{E,r,t)27rr- J dtj drfN,{E,r,t)27rr} 

poo poo rTo poo 

-{ dtj drfN^{E,r,t)27rr+ dtj drfN^{E,r,t)27rr} (B3) 
Jto Jo Jo Jrq 



c 



The integrations are carried out in the domain 1 and 2, where domain 1 is Tq < t < oo and 
< r < oo, and domain 2 is < t < Tq and Rq < r < oo, respectively. In the following, we 
discuss how to find the appropriate values of Rq and Tq for separation of distant and nearby 
source components. 

Taking the three dimensional solution of the formula (A2), the total flux J{E) is given 
by the formula (A3) as 

e;-7-(i+5)/2 

^<^) = s^*(Szwr3w^^i(^ - - 1/21. (B4) 

in the case of a power-law source spectrum, where B is the Beta function. Then the energy 
spectrum extends with a power-law of the form £'-')'-(i+'5)/2 -vvithout cutoff, as is already 
known in the literatures (e.g. Berezinskii et al. 1990). 
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The first term in formula (B3) (domain 1) is given by 

wliere i? is the incomplete Beta function. The fraction of the flux from domain 1 to the total 
flux, i.e. (B5)/(B4), should be small enough for the justiflcation on the separation of nearby 
sources from distant sources. Numerical value of the fraction (B5) / (B4) is less than a few 
percent when Tq = lO'Vr and E > 2TeV. 

The contribution from the domain 2 (the second term in formula (B3)) is given by 

To I7'-7 



o_ 



Replacing the variable t with 

1 



X — 



i-{i-bEty-^' 

the integration (B6) is reduced to be 

C ^-7-(l+<5)/2 



where 



with 



4.f^" (inm-6W'^ . d.(.-l)^x-»-^e-, (B7) 

l_(l_6BTo)l-* 

_ RlbE{l-6) _ 1 Rq ^ 
ADoE^ 2^rJ ' 

R, = {2DoE'/{bE{l - 



The physical meaning of i?a is that it represents almost the average distance traveled by 
diffusion of electrons with energy E. 

The upper limit of the integration in formula (B7) is obtained by putting the lower 
boundary of x to be 1.0, i.e. Tq = 1/bE. The integral of (B7) is then given by the Confluent 
Hyper-geometric functions, as 

7 - h , V^iF,[-0.5 - 0.5, -xo] 2^xriF,[-'-^, 1.5, -xp] 

r[i.5 + ^] ng] ^ ' 

or in another form of the Confluent Hyper-geometric function U, as 

,-.or[I^]C/[I^,0.5,xo]. (B9) 
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If the numerical value of (7 — 2 + 5) /(I — 5) is an integer, the integration in formula (B7) 
can be performed without imposing any conditions on Tq. For example, if we take 7 = 2.4 
and 5 = 0.3, we have (7 — 2 + 5)/(l — (5) — 1.0. Then the integration of (B7) is carried out, 
and we have 

g ^-3.05 

4^-^*^°(47rL'o(0.7)6)V2 



where F is the incomplete Gamma function. 

Here we define the ratio F of the flux in domain 2 to the total flux as 

where we used 5 [2, 1/2] = 4/3 for the total flux from formula (B4). F is a function of xq, 
and gives the relative contribution to the total flux for a given value of bETo from sources 
in domain 2. The numerical values of F is listed in Table 3. As shown in Table 3, we need 
to take Xo > 0.8 to reduce the contribution from sources in domain 2 to less than 6% of the 
total flux. This means we need to take the boundaries of domain 2 as i?o > V^Ra.- The 
average travel distance i?a is evaluated as 0.5 — O.Skpc at £^ = 2TcV depending on the values 
of Dq here we adopted. Thus if we put Ikpc for Rq, the contribution from domain 2 gives 
less than 6% for E > 2TeV. 

In summary, the treatment of separation of nearby and distant sources is justified in 
the energy region beyond ~2TeV if we put i?o=lkpc and Tq = lO^yr. 



Table 3 
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Fig. 1. — Energy loss coefficient of cosmic-ray electrons in the Galaxy with energy. The 
magnetic field is assumed to be B± = S/xG. 
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Fig. 2. — Contours of the electron flux at 3 TeV between distances and ages with the values 
of each SNR, in the case of the prompt release of electrons after the explosion. Lines show 
equal flux contour for Fq = (£'/GeV)^ J(GeV^m~^s'~^sr~^), where J is the flux of electrons at 
3 TeV. Contour levels are 10^ 10^, 10^ and 10° (GeV^m-^s-isr-^) with the output energy 
of electrons over IGeV of 1^ = 1 x 10^^ erg/SN. Here, the injection spectrum is a power- 
law with an exponential cut-off of £"0=10 TeV (dotted line), 20 TeV (dashed line), and 00 
(solid line). Upper and Lower panels show the flux contours for D — Do[E /TeV)^'^ with 
Dq = 2 X lO^^cm^s"^ and Dq = 5 x lO^^cm^s"^, respectively. 
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Fig. 3. — Calculated energy spectra of electrons without a cut-off of the injection spectrum 
for the prompt release after the explosion {t — 0), compared with presently available data. 
Here we took the diffusion coefficient of D = Dq{E /TeV)'^-^ with Dq = 2 x lO^^cm^/s and 
Do = 5 X lO^^cm^/s in TeV region, and D — 2xlO^^(£'/5GeV)°-^ in GeV region as given by 
the formula (3). See text for details. 
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Fig. 4. — Calculated energy spectra with a cut-off of -Bc^lO TeV and 20 TeV for the prompt 
release after the explosion (r = 0). Here we took the diffusion coefficient of Dq = 2 x 
lO^^cmVs. 
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Fig. 5. — Calculated energy spectra with D^ — Iy. lO^^cm^/s and a cut-off of Ec—20 TeV 
for the burst- like release at r = 5x10^ yr, 1x10^ yr, 5x10^ yr, and 1x10^ yr, respectively. 
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Fig. 6. — Calculated energy spectra with Dq = 2 x lO^^cm^/s and a cut-off of Ec=20 TeV 
for the continuous release time of r = 0^1x10^ yr and r = 0~lxlO^ yr, respectively. 
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Table 1. List of nearby SNRs. 



SNR 


Distance (kpc) 


Age(yr) 


£^max(TeV)- 


Reference 


SN185 


0.95 


1.8 X 10^ 


1.7 X 10^ 


(Strom 1994) 


S147 


0.80 


4.6 X 10=^ 


63 


(Braun et al. 1989) 


HB 21 


0.80 


1.9 X 10^ 


14 


(Tatematsu et al. 1990); 










(Leahy & Aschenbach 1996) 


G65.3+5.7 


0.80 


2.0 X 10^ 


13 


(Green 1988) 


Cygnus Loop 


0.44 


2.0 X 10^ 


13 


(Miyata et al. 1994); 










(Blair et al. 1999) 


Vela 


0.30 


1.1 X 10^ 


25 


(Caraveo et al. 2001) 


Monogem 


0.30 


8.6 X 10^ 


2.8 


(Plucinsky et al. 1996) 


Loopl 


0.17 


2.0 X 10^ 


1.2 


(Egger & Aschenbach 1995) 


Geminga 


0.4 


3.4 X 10^ 


0.67 


(Caraveo et al. 1996) 



^ Maximum energy limited by the propagation of electrons in the case of the prompt 
release after the explosion. The delay of the release time gives the larger value. 



Table 2. Ratios of the three dimensional solution without boundaries J3 to the two 
dimensional solution with boundaries J for various energies and halo thicknesses. 



Electron Energy 



/i(kpc) 


IGeV 


lOGeV 


lOOGeV 


ITeV 


lOTeV 


1 


3.35 


1.69 


1.08 


1.00 


1.00 


2 


1.77 


1.10 


1.00 


1.00 


1.00 


3 


1.30 


1.01 


1.00 


1.00 


1.00 


4 


1.12 


1.00 


1.00 


1.00 


1.00 


5 


1.05 


1.00 


1.00 


1.00 


1.00 
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Table 3. Relative contribution F of the flux at E = jhr in domain 2 to the total flux. 



0.2 


0.4 


0.6 


0.8 


1.0 1.2 


1.4 


1.6 


1.8 


2.0 


F 0.2722 


0.1481 


0.0901 


0.0580 


0.0387 0.0265 


0.0185 


0.0131 


0.0094 


0.0068 



^ <t <To, Ro <r <oo. 

^ See text for the definition of the domain and the parameter xq. 



